0-brane Quantum Chemistry 



VACLAV KARES^ 

Institute for Theoretical Physics, Masaryk University, 
- - - _ Kotldfskd 2, 611 37, Brno, Czech republic 

o 
o 

Abstract 

I ^ \ We apply two different numerical methods to solve for the boundstate 

^ I of two 0-branes in three dimensions. One method is developed by us in this 

■ work and we compare it to a method existing in the literature. In spite of 

considering only three dimensional Minkowski space we obtain interesting 
^ _ results which should give some basic understanding of the behaviour of 

Q\ • 0-branes. 



o 
o 



1 Introduction 



Oh! 

O . D-branes are very important non-perturbative objects in string theory. Their 

existence is essential since they are needed, among other things, for various string 
theory dualities to work. The most primitive definition of a D-brane is as a hyper 
^ ' surface on which open strings end. 

' In this work we are interested only in DO-branes. They play a particularly 

important role since they are the basic constituents of Matrix theory (Ij, which, 
being a suggestion for a non-perturbative definition of string theory (or M-theory 
[2]), deserves particular attention. It is also possible to think about them as 
bound states of higher dimensional unstable D-branes [3 E]. Furthermore, as 
was shown by Myers 0, they can also form bound states with all the properties 
of higher dimensional branes. It is therefore clear that DO-branes have many 
interesting properties that make them worth studying in detail. Concretely this 
means studying supersymmetric quantum mechanics which is an interesting topic 
in itself. Although this topic has been studied before we feel that there 

are still unresolved issues. In particular one could use a computer to compute 
properties of DO-branes which are not possible to address analytically because 
of the complexity of the theory. This is the goal of this work, to develop a "0- 
brane Quantum Chemistry". It should be mentioned that similar issues have 
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been addressed in jH]. However, we use a different method which we compare to 
the method developed in (to which we also suggest certain improvements). 

More concretely, in this work we try to find the bound state of two DO-branes 
in three dimensional Minkowski space (which is really a toy-model for the real 
situation, DO-branes in ten dimensional Minkowski space). It should also be 
noticed that the really interesting cases where the DO-brane theory is thought 
to describe macroscopic supergravity states is achieved by taking the number of 
zero branes to infinity. 



2 The model 

The low energy physics of parallel Dp-branes is governed by the dimensional 
reduction of 9 -|- 1 dimensional A/" = 1 supersymmetric Yang-Mills theory with 
U (N) gauge group to p+1 dimension ^H]- The center of mass motion is governed 
by the overall U{1) factor so if we are interested in relative motion only, we can 
choose the gauge group to be SU{N). In the case we are interested in, the relative 
motion of two 0-branes, we thus choose the gauge group to be SU{2). The action 

ism 

9s ^9s ^ 

and the Hamiltonian derived from this is [Zj 

^ ^9s j=i 

-l^^'-X^s {x'aXa - x'aXa) - \e'^''X^, {x\xa + X^a^a) (2) 
together with the constraint one gets from varying with respect to Aq 

^ e'^'^X^T^t - itAX'^A) = ^ ■ (3) 
On quantum level, we restcrict our Hilbert space to vectors which satisfy 

G"|^) = (4) 

that is, our physical space is gauge invariant because are gauge generators. 

We will only study motion of two 0-branes in three dimensional Minkowski 
space. We hope that this gives us the basic behavior of 0-branes and also an 
understanding of the full problem. This problem of two branes is also described 
in [7] but we study it in a different way. In the three dimensional case the action 



2 



is given by dimensional reduction of A/" = 1 supersymmetric Yang-Mills theory 
with SU{2) gauge group in 2 + 1 dimension to + 1 dimension. The Hamiltonian 
takes the slightly simpler form 

-^e'^'^X^, {xY - X'X j - ^e^'^^X^ {xY + x't) (5) 

where Xj are fields in the SU (2) adjoint representation and x is a complex fermion 
also in the adjoint representation. Of course we still have to impose gauge in- 
variance (jlj. In fact, the gauge invariance complicates things somewhat since we 
would like to separate out gauge invariant degrees of freedom from pure gauge 
degrees of freedom in our basic quantum mechanical operators and tt°'. Let 
us focus on physical content of the X". It contains six components (the gauge 
index runs over three values and the space index i runs from 1 to 2). We know 
that we can remove three of these variables using gauge transformations so only 
three variables are observable. These three variables should describe the relative 
position of two pointlike objects in two space dimensions. We draw the conclusion 
that one of the physical variables do not have the interpretation of a coordinate 
but rather as some internal auxilliary degree of freedom. 

To get some further insight into this problem it is neccesery to investigate 
the bosonic vacuum of the theory. It is possible to explicitly separate the gauge 
degrees of freedom from X^ " by decomposition it in matrix form 

(X),, = {^)ar{A)rs{v)si • (6) 

Here the matrix ip is an group element in the adjoint representation of SU{2). 
Thus when the gauge group acts on Xf , ip just changes by ordinary gauge group 
multiplication (from the left). We will parametrise the group element ip by the 
"angles" 7 ()A.5|) . This decomposition has the advantage that all the gauge 
dependence sits in ip and all the other matrices are gauge invariant. 

In an analogous way we have separated out the dependence on rotations in 
space. Namely, performing an SO {2) rotation in space we have an element of 
5*0(2) acting from the right on the matrix Xai- Thus we can separate out the 
dependence on the angle in space (we will call it 0) by saying that r/ is a group 
element of SO (2). 

We are left with the matrix A ()A.2|) which by construction is both gauge and 
space rotation invariant 

/Ai \ 

A = A2 . (7) 

Vo 0; 

The bosonic potential in is gauge and rotation invariant and in the new 
decomposition coordinates depends only on two Aj which have length dimension 

(Fig. m 
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Figure 1: Bosonic potential 



The parametrisation 



Ai = r cos 6; X2 = r sin 6 



(8) 



is the only way how to obtain exactly one variable, r, with the dimension length, 
which could represent relative distance of two branes. The dimensionless 6 is the 
auxilliary coordinate. The potential in this coordinate reads 



Looking at the picture we can draw some interesting conclusions. If we fix 
a point on the bosonic vacuum (a classical static configuration with minimum 
energy), that is on the axes, we can study the behavior of the potential for small 
fluctuations of the auxilliary variable 6. We see that for large r the 6 fluctuation 
are very much suppressed but at small r, 9 will be allowed to fluctuate. This can 
be interpreted to mean that when the branes get close to each other, they can 
start to move also in the 6 direction. Thus, 6 is an auxilliary coordinate which is 
visible only when the branes come close together. 

The above discussion included only the bosonic degrees of freedom, we should 
keep in mind that the fermionic degrees of freedom can (and will) change this 
behavior somewhat. In essence, the Pauli repulsion will try to spread out the 
wavefunction as much as possible. 

When we use the standard operator representation of vrf and the first term 
in the Hamiltonian (jSJ is proportional to the Laplacian which we have to rewrite 
in the decomposition coordinates above to separate out gauge using formula 



where the g is the metric which is trivial in coordinates. It is also good idea 



— sin^ 26 . 

8gs 



(9) 
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(10) 
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to rewrite the Laplacian in terms of gauge angular momenta L° = e'^^'^X^irf: 



^1 d . . d . d 

L = — zcot7smQ;7r — h «csc7sma— + zcosa- 
aa op OJ 

t2 ■ d . d . . d 

L = cot 7 cos 07- — h «csc7cosa— — zsma— 

Oa op O'j 

d_ 

da 



and physical angular momentum 



(9 



(12) 



since their action on states with given gauge and rotational properties is simple. 
In particular, the (bosonic) ground state should have total spin equal to zero and 
be gauge invariant. All angular momenta are Killing vectors of the metric. 

In Appendix^ we derive the lagrangian expressed in terms of gauge invariant 
variables and angular momenta 



L 



-r^— + ^^— sin4^— 



l,T.-l^ 



(13) 



Here L^^ fi = 0, 1, 2, 3 is a compact notation for the (physical and gauge) angular 
momenta defined above. Furthermore, we have defined 



-1 



(14) 



and n 



sin 6 



cos^ 6 



\r^sin26' 



' sin 26 \ 



J 



(15) 



To find the (bosonic) ground state we will find all gauge invariant states with 
spin zero. The first state is the vacuum state |0). Then we may act with the 
fermionic creation operators on the vacuum to find new states. The following 
states has total spin zero and they are gauge invariant 



|r) = i^ar |a) = V^A.^e'^VxlO) . 



(16) 



That is, they satisfy 



|r) = . 



(17) 
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The most general gauge invariant wavefunction with total spin zero can then be 
written 



g{r,e)\Q) + Ur,e)\r) . (18) 

It would also possible to construct the superpartner ground state with the help 
of 

\r') = e'HarX^ |0) 

|0')= 1 2^WxV|o) . (19) 

Now we study what happens when we act with the Laplacian pHj) on this 
wavefunction and using the result we write the Hamiltonian matrix in this base. 
Let us start with second part of the Laplacian which is the relevant bosonic piece 

- (M-i*-i)„bM^ |r) = \s) - |r) TrR-^ 
-2(M-^*-i)oaLL" |r) = 2i|M)no>„.. 
-(M-i*-i)ooLL|r) = -|r)noV (20) 

where Tr is only on the gauge indexes. Finally the fermionic interaction in the 
Hamiltonian gives us 

Ep\r) = -|0)(Ar7),i + i|0)(Ar7),2 

Hf\0) = -\r){Ar]U-i\r){Ar])r2. (21) 
So the full Hamiltonian matrix elements are 

Hor = -(Ar^),! + i(Ar/)r2 

^^00 = -^h + ^rUm^29 (22) 
where 

Notice that the Hamiltonian we have obtained is the same as in [Hj but differs 
from the one used in [7j. 



3 Numerical calculation I 

This section gives an overview of solving the Hamiltonian eigenvalue problem (0) 
using the numerical renormalized Numerov method [S]- With it one can solve 
for the discrete spectra of a one dimensional operator of the form 

-i^i + vw (24) 
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which acts on L[a, 6] with Dirichlet boundary conditions. However, our 
Hamiltonian depends on the two coordinates r and 6 (the other coordinates, the 
angles, are fixed by the requirement that we are studying only gauge invariant 
states with spin zero). The "kinetic" term in our Hamiltonian (i.e. the term 
which contains the derivatives) is 

or or r^ sm 46 oa Oo 

and we see that it is naively not of the form required above. Let us sketch 
briefiy how to modify our problem to be able to apply the method. An arbitrary 
wavefunction can be written in the form 

l^) = E%(0>^..(^)e. (26) 

where {Yij{9),j = 1, . . .} is a complete basis of functions (with apropriate bound- 
ary conditions) in 6 (the explicit choice of basis does not have to be the same 
for different values of the index i but can be chosen to optimize the numerics). 
The functions {"^ijjj = 1, . . .} (which we will compute by the Numerov method) 
one can think of as being "combination coefficients" depending continuously on 
r. {cj} is the standard base spanning C", in our case it is the four dimensional 
complex vector space on which the Hamiltonian matrix acts. Choosing a con- 
crete basis depends only on the Hamiltonian domain which depends on one of 
the self conjugated extensions of the Hilbert space (an extension of the Hilbert 
space such that the Hamiltonian operator is hermitian) but we will rather apply 
a physical principle which will be described later. The expectation values of the 
Hamiltonian in the base 

YijCi (no sum) (27) 

gives the same number of coupled equations for the radial part as the number 
of in ()26|) . Thus the problem is now correctly defined and the matrix repre- 
sentation of the expectation value of the Hamiltonian in the basis above forms 
the potential which is used in the one dimensional Hamiltonian (j24|) . To get 
1 in the kinetic term we have to orthonormalize the base. Furthermore, to be 
able to use the Numerov method on a computer we need a finite basis which 
means that we need to "cut off" or restrict the base to be finite. The rescaling 
of the wavefunction |\E') by the factor r^/^ transforms the radial part of (|25|) to 
the operator 

and changes boundary condition to Dirichlet. 
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Let us now describe the physical principle we use to choose the basis functions 
Yij{6) for our Hamiltonian (0). 



We define the wavefunction on the interval 9 G [0, 7r/4] (jA.27|l . It is not 
necessary to choose this particular interval, one could, for instance, select the 
interval [— 7r/4, 0] instead of the above mentioned. Using the identification 6 = —9 
the Hamiltonian defined on the interval 9 G [— vr / 4, 0] acting on states with total 
spin zero is connected to our original Hamiltonian on the interval 9 G [0, 7r/4] by 
the unitary transformation 



H{9) = U^H{9)U 



(29) 



where 



U 



/I \ 

1 

-1 

\ 1/ 



(30) 



It is also possible to consider other interval 9 G [7r/4, 7r/2]. Using the identification 
9 = 71/2 — 9 the corresponding Hamiltonian can be also obtained by a unitary 
transformation with the matrix 



U 



/I 



V 



1/ 



(31) 



The wavefunctions of course also transform under the unitary transformation 



\^ili9)) = U\^i9)) . 



(32) 



If we require that the wavefunctions be everywhere smooth, the above condition 
severely restricts the possible wavefunctions and in particular the basis wavefunc- 
tions that we can use. We note that this principle was not used in [7] and hence 
the 9 derivative of their groundstate wavefunction at ^ = 7r/4 is not well defined. 
For the bosonic part of the wavefunction we will use the basis states 



cos(4j6')ei j = 0, . . . 
and for the fermionic part we will use 

- i cos 6* cos(2j6')e2 + (-1)^' sin 9 cos(2j6')e3 j = 0, 



(33) 



(34) 



They nicely cancel the divergences in the Hamiltonian coming from the Laplacian. 
Of course these basis functions satisfy the boundary condition above and they 
also form a complete basis for the functions with the given boundary conditions. 
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However, the basis is not orthonormal due to the non trivial 6 part in the measure 
()A.26|) . So we need to orthonormahze these basis functions to get 1 in ()24|) . 
Doing so, using the Gramm-Schmidt procedure, we diagonahze the 6 part of the 
Laplacian. Notice that these two sets of basis functions also correspond to the 
bases which are implicitly used in the Wosiek method [Hj when one are doing 
calculations with states with total spin zero (the basis functions above depend on 
the total spin). This will be shown in sectional Thus we will be able to compare 
our results and our method very directly with the results obtained by the Wosiek 
method 0. 

It is not obvious how much our results for a fixed number of basis functions 
fit the exact solutions which one would get using the complete basis. To get 
some intuition for how the general solution would look like we will repeat the 
calculations increasing the number of basis functions each time and hopefully 
one can extrapolate the result to the exact case. At least we should be able to 
make an intelligent guess at the properties of the exact solution. 

To study the groundstate of our Hamiltonian (which, because of supersym- 
metry, should have energy zero JT]) is a good test for the method described 
above. 

There are some results for coupling constant gs = 0.1 in the table [T] (p. ITTjl 
where N is the number of the test functions ()33p34j) and E is their corresponding 
groundstate energy. The dependence of the energy on the number of functions 
we include in the basis given in the last two columns of the table is as follows 

1.44 , , 

E = ^. (35) 

As was claimed in jB], one can therefore predict the groundstate energy as a 
function of the number of basis functions with very high accuracy. We therefore 
see that in any concrete numerical calculation (using a finite basis) we do not 
expect to get zero energy. Only in the (numerically unobtainable) case of infinite 
basis do we get zero energy. 

The picture (Fig. ^ is the probability density for the case with the high- 
est number of basis functions in the table. The following two pictures show the 
bosonic contribution coming from ()33|) and the fermionic one coming from ()34|) . 
The domain of these plots is (r, 6) G [0, 2.2] x [— 7r/4, vr/4]. The hill of the prob- 
ability density is located at the boson potential valley Q and isolines represents 
sections for fixed r, 6 and one for fixed density on each picture. The maximum of 
the probability density of any constant r section is in the potential valley (^ = 0). 
Notice that the global maximum is not at r = 0. This is probably an effect of the 
fermion Pauli repulsion which can be seen from the purely fermionic contribution 
to the probability density which is zero at r = 0. 

Increasing the number of basis functions, the only thing that happens is that 
the global maximum moves slowly to larger and larger r at the same time as the 
whole wavefunction becomes more spread out in r but more peaked in 6. One 
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Figure 2: Probability density 




Figure 3: Bosonic part of the probability density 
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Figure 4: Fermionic part of the probability density 



can expect that considering the complete base ^T7\i the ground state density will 
be the same near the origin as in the picture and also will have the hills on the 
valleys of the potential which will be sharper and sharper when we follow the 
potential valley to large r. One can not conjecture in this case the behavior of 
the hills by this numerical method. Rather one has to use other methods J2] for 
the asymptotic behavior of the wavefunction at large r. 



4 Numerical calculation II 

Here we briefly review the body of paper jB] where a different approach to solving 
the eigenvalue problem of a Hamiltonian is given. The main idea is to choose a 
finite subset of a complete basis of the Hamiltonian Hilbert space and then find 
the combination which minimizes the lowest lying energy for instance. A lot of 
Hamiltonians have the form (|24j) which can be easily rewritten in terms of the 
harmonic oscillator 

-JS^i + jCV (36) 

creation and annihilation operators a+,a_. This is true even in the case where 
the potential V{x) is not quadratic. The complete base above is formed by eigen- 
vectors of this linear harmonic oscillator which are constructed by creation oper- 
ators acting on the vacuum. If we choose a finite subset of this basis, calculating 
Hamiltonian matrix elements is very easy. 

In the harmonic oscillator Hamiltonian which we used to define the ba- 
sic creation and annihilation operators there is an arbitrary parameter ( in the 
potential which we may use to optimize the numerics. The option to use ( to 
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optimize the calculations was invented by us and should be regarded as sugges- 
tion for an improvement of the method presented in jH]. Concretely we will do 
it like this: first we choose a subset of the basis functions, we then compute the 
energy (or whatever else we would like to calculate) for different values of (. We 
then find the value of ( which gives the minimum value for the energy. It turns 
out that it is much more effective to minimize the energy in terms of ( than to 
try to use a larger basis since the computer time used for a calculation increases 
exponentially with the number of basis functions included. 

There is one extra complication in our problem. Because of the gauge sym- 
metry one should also restrict the basis to include only gauge invariant states. 
The generalization to the fermionic part and to the higher dimensional problem 
is straightforward. 

Let us illustrate the method above on our Hamiltonian. First of all we have 
to form all gauge invariant states which can be obtained by acting with a com- 
bination of creation operators on the vacuum. The SU{2) algebra has only two 
tensors which can form scalars from vectors 



These vectors are formed by creation operators a+^, x"" and a combination of 
them has to be contracted in gauge indexes with the special tensors above to get 
a gauge scalar. Let us write all non zero fundamental possibilities 



where any combination of them acting on a gauge invariant state gives also a 
gauge invariant state. One has to start with the vacuum to generate all gauge 
invariant state. However, if we are looking for the groundstate we are interested 
only in states which can form a state with total spin zero. That is only states 
which consist of the fermion vacuum or two creation fermion operators acting on 
it because these are the only cases where the fermionic contribution to the total 
angular momentum can be canceled by the bosonic angular momentum. So using 
the operators 



6' 



(37) 




(38) 




(39) 



(40) 
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gives the appropriate complete basis for our Hamiltonian. Choosing a finite subset 
of the basis states allows us to calculate matrix elements of any operator on the 
computer. It is a good idea to orthonormalize the states in the basis to avoid 
problems in the eigenvalue calculation when the metric is non trivial. See the 
next section for the explicit form of first few states in the basis. 

We are able to illustrate the effect of the optimized ( now. Let us show what 
happens with the ground state energy when we calculate with fixed ( but with 
different coupling constants Qs G (0, 2] in our Hamiltonian (0). In this calculation 
we have chosen the restriction of the basis to include up to six creation operators 
in the bosonic basis and up to five creation operators in the fermion part. The 
result is given by the dotted line in the figure El 




Figure 5: Groundstate energy dependence on the gs 



It seems that the energy of the ground state diverges when the coupling con- 
stant goes to zero which is clearly wrong! The way out of this dilemma is to 
remember that we also have the parameter ( at our disposal. If we for each gs 
use the ( which gives the minimal energy, we get the dependence of the ground 
state energy on the coupling constant as drawn in bold line. This is clearly con- 
sistent with the g^^ dependence one expects from dimensional arguments. This 
is a very nice illustration of the power of our suggestion for improving the Wosiek 
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method jHj by including the ( parameter. One can guess that there has to ex- 
ist a couphng constant such that the optimahzation is not needed. This is also 
viewable from the picture. 

There are also some results for coupling constant Qs = 0.1 in the (Tab. ^ p. 
fTTj) where the meaning of the columns is described above the table. The energy 
of the bold line for = 0.1 corresponds to the first two columns in the table 
for the highest order obtained by the Wosiek method. We draw the density plot 
(Fig. ini) for this calculation on which we can see the bump which for a larger 
number of basis functions develops into the long and thin hill in the picture on 
pageimi Its domain is {r,6) e [0,1.5] x [-7r/4, 7r/4]. 




Figure 6: Groundstate by Wosiek method 



5 Comparison 

To compare our method with the Wosiek method it is important to know the 
correspondence between the base (j33l34|) and the base formed by creation oper- 
ators. Since we have written our Hamiltonian in terms of angular momentum 
operators it is useful to classify all states in terms of angular momenta. In fact, 
this is even more useful since when we are discussing the ground state with total 
spin zero. We may rewrite the angular momentum operator (|T^ as 

i(a4a4-a+>4) (41) 

using the correspondence between creation/annihilation and Xf,7rf operators 

a±t = ^{X: T i/C<) • (42) 
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Let us write down how the lowest lying angular momentum states look like. First 
we do the purely bosonic sector. The (zero angular momentum) ground state |0) 
in the coordinate representation looks like 

3 

V'o(r) - e-^f^l (43) 

Next we can write down all gauge invariant states created by two bosonic oscil- 
lators. 

|1') = -^a+W?|0) 
|2') = i=a+X2|0) 

|3') = ;^«+>+2 |0) . (44) 
The following linear combinations 

i 11') -5 12') + ;^ 13') 

4 11') + ^ 12') + ^13') (45) 

are orthonormal eigenvectors of the angular momentum operator with eigenvalues 
{—2,0,2}. Using formulas for the creation/annihilation operators in the X",7r" 
representation we may write the coordinate representation of these states as 

^r\cos2^e-^^Vo 

i=(-3 + r\)V'o 

— ^r^Ccos2^e2^% . (46) 

The angular momenta spectrum of the next six states all formed by four creation 
operators is {—4, —2,0,2,4} © {0} where the two orthonormal states with zero 
angular momenta are given by 

^ (60 - 40r^C + + 3r^C^ cos 4^)V'o 



4a/105 



^ (36 - 24rX + 3r\2 _ ^4^2 ^q^^^ _ ^47^, 
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They give the first non trivial wave function dependence on 6 in the zero angular 
momentum sector. The 6 part in higher order states with zero angular momenta 
reproduce the basis (j33|l which we guessed as the Fourier basis for the functions 
with the correct boundary conditions. We see that a new cos 4j6 appears when 
there appears a new in the angular momenta spectrum. 

For reference we also give the spectra for the states created by six and eight 
creation operators {-6, -4, -2, 0, 2, 4, 6}©{-2, 0, 2} and {-8, -6, -4, -2, 0, 2, 4, 
6, 8} © {—4, —2, 0, 2, 4} © {0}. With this we close the discussion about the purely 
bosonic part. 

In the same way we can study the fermionic part. The wave function that 
interest us most are the ones with angular momentum one since they are the 
ones which contribute to the ground state with total angular momentum zero. 
On the lowest level in the creation/annihilation operator basis we have two gauge 
invariant orthonormal states 



which are eigenstates of angular momenta with eigenvalues { — 1,1}. It is again 
possible to rewrite these states in the coordinate basis with help of the orthonor- 
mal and gauge invariant purely fermionic states (|16|) as 



The 6 part of the second state gives the first function in the basis ()H4|1 . The 
states in the next level consisting of the six states we get by acting with all gauge 
invariant combinations of two bosonic creation operators on the previous two 
states have angular momentum spectrum {—3, —1, 1, 3} © {—1, 1} where the two 
orthonormal states with angular momentum one are given by 





They can be combined into the orthonormal states 




(49) 




(50) 
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(51) 



This agrees with the second function in our fermionic basis. Continuing we find 
that the 9 part of the angular momentum one states at higher levels reproduce the 
basis ()34|) which we have used. Writing the angular momentum spectrum for the 
next level (states obtained by acting with all gauge invariant combinations of four 
bosonic oscillators on the original fermionic states) we have {—5, —3, —1, 1, 3, 5}© 
{—3, —1,1, 3} ©{ — 1,1}. This closes the discussion about the fermionic part and 
again gives a hint how the general spectrum looks like. 

We have now achieved our main goal to rewrite the states in the Wosiek 
method in terms of angular momentum eigenstates in order to be able to compare 
the two methods. Since we now know at what level a particular basis function 
appears in the Wosiek method, we know where we should "cut off" the set of 
basis functions used in our method in order to get an equivalent calculation. 
Furthermore, doing the calculation using the two different methods at the same 
level we should still expect that our method should give better results since in 
our method the radial part of the wave function is not a priori fixed while in 
the Wosiek method the radial part comes together with the angular part, see for 
instance ()47|) . Also, our method should be more effective for fixed spin calculation 
on the computer since we have the possibility to use only the angular momentum 
eigenstates we need while in the Wosiek method one is using all possible states 
which leads to a huge number of states not contributing to the wavefunction but 
absorbing space and time in the computer. 

To illustrate this we compose the following table. We compare the (improved) 





c 


Ep 


N 


E 


4+2 


4.55 


0.395 


1 + 1 


0.376 


10+8 


4.55 


0.303 


2+2 


0.295 


20+20 


4.73 


0.216 


2+3 


0.214 


1540+440 






10+10 


0.075 


11480+3080 






20+20 


0.037 


37820+9920 






30+30 


0.025 








40+40 


0.018 








50+50 


0.014 








80+80 


0.008 



Table 1: Results 

Wosiek method with Np boson and fermion creation/annihilation states of the 
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type ()44|48|) with our method using only the test functions ()33I34|) to the cor- 
responding level for = 0.1. The ( given in the table is the optimized value 
for the (improved) Wosiek method, the N counts the number of the test func- 
tions corresponding to the level used in the Wosiek method and the groundstate 
energies obtained by both methods are also given. 

We see that indeed the energy calculated using our method is lower than in 
the Wosiek method. Furthermore it is clear that the number of states used in 
the Wosiek method increases uncontrollably as compared to the increase of basis 
functions used in our method. 

Finally we also compare probability density sections for ^ = in the case of 
the third line in the table (the last case where both methods have been used) 
on the domain r G [0,2]. The states which generate these sections have the 
same norm to compare them correctly. We have also included higher order wave 
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Figure 7: Comparison of wavefuctions at 6' = 

functions calculated using our method. The labels indicate which order they 
correspond to. These wave functions are not equally normalized but have rather 
been chosen to have the same value at r = in order to illustrate the behavior of 
the wavefunction when the order is increased. We see that the maximum moves 
outwards at the same time as the wavefunction becomes wider. 
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6 Conclusions 



The aim of the present work was to find an approach how to calculate the bound- 
state of two 0-branes. This we have achieved with good results. In the course 
of our work we also found a way to improve the Wosiek method but our method 
still gives better results. As this work was being completed the closely related 
paper J7j by Wosiek (which improve his original method) appeared. 

We have also calculated the groundstate probability density near the origin 
with high accuracy which gives a basic intuition about the physics of branes on 
the string scale. In particular about the meaning of the auxilliary coordinates 
which become important at small distances between the branes. It is for example 
interesting to observe that the most probable position of the branes is not on top 
of each other but rather at some small distance away from each other. This we 
understand as an effect of the Pauli fermionic repulsion. 

It is possible to apply the method described here for branes in higher dimen- 
sional Minkowski space. However there are a lot of complications which do not 
allow numerical results of high precision. Namely, there are many more auxilliary 
coordinates in the higher dimensional case (or in the case where there are more 
than two branes). Since our basis functions are functions of the auxilliary coor- 
dinates, this implies that our basis functions will be multi dimensional functions 
and thus the number of components that the computer has to handle increases. 
Naively speaking the number of effective radial wavefunction in the state com- 
puted by the Numerov method for a given dimension is the number in our toy 
model raised to the power of the number of auxilliary coordinates. This fact of 
course taxes the computer very heavily. Another important point that we wanted 
to study was how the branes behave in external (generalized electric-magnetic) 
fields. The potential coming from these background fields represents this inter- 
action. However, our two dimensional toy model is too simple to display any 
interesting effects. Again one would have to go to higher dimensional Minkowski 
space to do be able to do some interesting calculations. 

As a first attempt to go further it would therefore be nice to consider three 
dimensional Minkowski space which would allow more general background fields 
which coming from other types of branes. However, we do feel that doing realistic 
calculations of the full ten dimensional (or even four dimensional) theory or in- 
cluding more than two 0-branes would not be possible even on a super computer. 
This is one of the reasons why we have not pursued this issue further. 
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A Metric derivation 

It is possible to decompose |12] in matrix form as © 

X = tpAr] (A.l) 

where ip is a. group element in the adjoint representation of the gauge group and 
?7 is a group element of the rotation group 5*0(2). The matrix A is then gauge 
and rotation invariant. We may now choose our physical coordinates as 




A = I A2 I , Ai = rcos^; A2 = rsin^ (A.2) 
together with 

r/ = exp(0T) , (A.3) 
where T is the generator of spatial rotation 



1 
-1 



(A.4) 



The matrix ip which describes the gauge degrees of freedom is parametrized by 
Euler angles a, /?, 7 

ip = exp (aTs) exp (7T1) exp (pT^) (A. 5) 

where the matrices are generators of SU{2) in the adjoint representation. 

The flat metric in the original coordinates is trivial and can be written in 
matrix notation as 

g = Ti{dX^dX) . (A.6) 
The decomposition (jH)) leads to the expression 

dX = di)K7] + %l)dKr] + %l)Mr] (A. 7) 

and then the metric is 

Tr(c/A^rfA - i>'^d^kk^i)'^di) + drf Mr] + 27]drf k^^'^ di)k) . (A.8) 

We see here that the differentials in dA are not mixed with the other ones and 
thus the metric will be block diagonal 



9i 

9m 

20 



(A.9) 



where the gi block is associated with the dK differential and Qm is associated with 
the differentials drj, dip. Now we compute the metric term by term. The A part 
in the case of two branes is 



Tr(dA'^c/A) = dXidXi 
or explicitly in variables r, 9 from ()A.2jl 

dXidXi = dr^ + r^dO^ . 



(A.IO) 



(A.ll^ 



Next we take care of the terms which contribute to the Qm part of the metric. 
Using the fact that 



the purely gauge part 



Tr(V'' #AA' i)' dip) 



(A.12) 



(A.13) 



can be written with help of the left-invariant one forms defined by the ip being a 
group element in the adjoint representation of SU (2) 

CUl = = C^LaT, (A. 14) 

where are the generators. In our coordinates we explicitly have the expressions 



uli = sin 7 sin pda + cos pd'j 
ujl2 = — sin J cos (3 da + sin pd'j 
i^L3 = cos 'yda + dp. 

Introducing the matrix Ql of the coefficient of the u 1-forms 

sin 7 sin /3 cos /3 
— sin 7 cos /3 sin j3 
cos 7 10 



(A.15) 





(A.16) 



n, 



the gauge part (|Xl3| IS 



(A.17) 



where = {da, df3, d'j) and the diagonal matrix Hat = Tr(TaAA-^Tf,) has ele- 
ments 



(A.18) 
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The calculation of the contribution from terms which contain drj is very simi- 
lar. Introducing the right-invariant one form defined by drji]'^ and using our 
parametrisation ()A.3j) we have 

Ti{drf Mr]) = r^d(j)'^ (A. 19) 

and 

TY{2ridr]'^A'^ij'^dijA) = 2r^ sin 29d(j)ujL3 ■ (A.20) 

In this expression it is necessary to write the ujl ()A.15|1 explicitly. Since we know 
that the metric is block diagonal ()A.9jl we will write the Qm part of the metric 
compactly. This is done by writing the matrix pSj) 



n 



/ r^sin26'\ 
sin^ 9 

cos^ 9 

\r'^sm29 / 



(A.21) 



^L={^ n ) • (A.22) 



and the structure matrix 

Then it is possible to write the sum of ()A.17IA.19l IA.20|) as 

e^nlnflie (A.23) 

where e-^ = {dcp, da, dj3, dj). So 

gm = nlun (A.24) 

and the full metric is 

dr^ + r^d9^ + e^fllUflLe . (A.25) 
It is now easy to find the square root of the metric determinant 

^r^ sin 7 sin 461 (A.26) 

which also gives us the integration measure. On the domain of (r, 9, 0, a, (3, 7) 

(0, 00) X [0, 7r/4] X (0, 27i] x (0, 2n] x (0, 2n] x (0, tt] (A.27) 

the transformation rule (jHI) is one to one and on its inside the Jacobian is nonzero 
and thus regular which is important in calculation of derivatives with respect to 
X, in terms of these coordinates. 
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B Laplacian derivation 



The metric in the new coordinates is ()A.25j) 

dr"^ + r'^de^ + e^rj^Hfiie (B.l) 

where e'^ = {d(f),da,dp,d'y), fli is defined by flA.22j) and 11 by (fT5|) . From the 
metric or directly from ()A.8|1 one can see that the metric has the block diagonal 
structure ()A.9|1 



9i{r) 



(B.2) 



where Qi is a diagonal matrix being the metric in the r, 6 variables and Qm is the 
metric in the other variables. Thus we are able to calculate the inverse of the full 
metric in terms of the inverses of gi and gm separately. Since gi does not depend 
on the gauge variables, the Laplacian also splits into two pieces 

1 „ . n ^ . 1 



, . , dA^^detg.detg^gra,) + -^^d,{^ det g,r.9^^d,) (B.3) 
Vdet det " V det gm " 

where the latin index runs over the {r, 6} variables and the greek indexes over 
the angular variables. Using the square root of the metric determinant ()A.26|) 
the first part of the Laplacian generated by gi is 

-j==d,{^deigg^d,) = -t^^V + 2 ■ ^zj ^^^^^^T^ • 

At this point we would like to rewrite the second part of the Laplacian in terms 
of the angular momenta ()12lll|) 

^, d d d 

L = — «cot7smQ;— h «csc7sma7— + zcosa— 

oa op Oj 

t2 ■ d . d . . d 

L = cot 7 cos a— — h «csc7COSq;— — zsma— 

Oa Op Oj 

. 4. ,B.5) 

We notice that the angular momenta are equal to the vector field Xr being dual 
to the right invariant one form ur = dipip'^ 

■ r .9 . d d 

— iLi = Xri = — cot7smQ;7- — h esc 7 sm a— + cos a— 

oa op c/7 

•r 9 d . d 

—iL2 = Xr2 = —cot 7 COS 07: — h CSC7 cosaT"— — sma — 

Oa op O'j 

-tL, = Xn3 = (B.6) 
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Since the right invariant one forms induce left translations, they are Killing fields 
of the metric. Also, since we know the relation between the left and right invariant 
one forms 

UR = #^-1 = ^LJL^-^ (B.7) 

we may transform the expression for the metric, which is expressed in terms of 
Xl into an expression in terms of Xr which is then equivalent to an expression in 
terms of the angular momenta. Explicitly we write the vector field Xl,r in terms 
of coefficient matrices Ql,r inverse to the matrices flL,R introduced in ()A.16|) . 
The relation ()B.7|) then implies 

Ql = ^P-'Qr . (B.8) 

To treat all the angles on equal footing, we incorporate also the physical angle 
into this formalism by extending the three by three matrices Qr, Qr, ip to four 
by four matrices by putting —1 in the "zeroth" row and column. These new 
matrices we denote by bold letters, for instance 



Now all vector fields obtained from are equal to — i times the physical and 
the gauge angular momenta and they are all Killing vectors of the metric. We 
can now write JIl = '^'""'^^^r; ©l = '^'""'^©r- 

We may now rewrite the angular part of the Laplacian in terms of the angular 
momenta which are given by ©r. The appropriate part in the Laplacian is 

= ©^M-i*-^©^ . (B.IO) 
First we expand the formula for the Laplacian 

-^d.i^g^'^d.) = ^d,{^)g^^d, + d,{gnd. + r^.Q. ■ (B.ll) 

The first term is reduced to 

^ :©«)/(M-^*-i©«)/^ (B.12) 



dx^' dxP 

where we have used the Killing equation for @r. In the same way the second 
term becomes 
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and finally the third part is 



@,)^^(m-'^-\,i&nV^^^ . (B.14) 



So the sum of these three contributions gives 



d 



:b.i5) 



where we have used that the angular momenta are defined in terms of the matrix 
Then the whole Laplacian is 

-^^'^ + -Y^A^'^^(^i-n - (*n-i*-i),.L''L^ . (B.16) 
or or r^ sm 46' o9 o9 
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